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Using a model Hamiltonian approach including electron-electron interactions, we systematically 
investigate the electronic structure and magnetic properties of chiral graphene nanoribbons. We 
show that the presence of magnetic edge states is an intrinsic feature of smooth graphene nanoribbons 
with chiral edges, and discover a number of structure-property relations. Specifically, we study the 
dependence of magnetic moments and edge-state energy splittings on the nanoribbon width and 
chiral angle as well as the role of environmental screening effects. Our results address a recent 
experimental observation of signatures of magnetic ordering in chiral graphene nanoribbons and 
provide an avenue towards tuning their properties via the structural and environmental degrees of 
freedom. 
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Graphene and derived nanostructures exhibit a large 
number of novel electronic properties P, Q • One of such 
features is the presence of electronic states localized at 
the edges of this two-dimensional (2D) nanomaterial ■s'] . 
These zero-energy edge states were predicted to give rise 
to a novel type of magnetic ordering Q which may lead 
to practical carbon-based magnetic materials [Hj. Even 
more appealing is the prospective of realizing the the- 
oretical proposals of novel spintronic devices based on 
graphene [6l4lCll|. While evidence for edge states has 



been seen experimentally [ll|, [12| , no solid proof of low- 



dimensional edge magnetism in graphene was presented 
until now. 

A recent scanning tunneling microscopy /spectroscopy 
study of graphene nanoribbons (GNRs) with ultrasmooth 
edges showed the presence of edge states with characteris- 
tic splitting in the dl /dV spectra - an unambiguous indi- 
cation of magnetic ordering [l^. These GNRs produced 
by the chemical unzipping of carbon nanotubes are 
chiral, i. e. characterized by low-symmetry orientation of 
the edges rather than by high-symmetry zigzag and arm- 
chair directions. While the presence of edge states at the 
chiral graphene edges is broadly recognized 
theoretical investigations of magnetic ordering driven by 
electron-electron (e-e) interactions have so far focused 
only on high-symmetry zigzag edges. 

In this Letter, we systematically study the electronic 
structure of chiral GNRs using a self-consistent model 
Hamiltonian approach including e-e interactions. In par- 
ticular, we investigate GNRs characterized by a broad 
range of chiralities and widths as well as address the 
effects of varying e-e interaction strength. Our study 
reveals that spin-polarized edge states are an intrinsic 
feature of chiral GNRs, in agreement with the recent ex- 
perimental observations. Moreover, we find a number 
of structure-property relations and unambiguous signa- 
tures of magnetic ordering of edge states which open new 



prospectives for their further exploration and for devel- 
oping practical spintronic devices based on them. 

The method employed in our study is based on the 
mean-field approximation to the Hubbard Hamiltonian 

+ {ni^{nii) + {nif)nii - {ni^){nii)) . (1) 

i 

The first term is the nearest-neighbor tight-binding 
Hamiltonian in which Cia (cl^) annihilates (creates) an 
electron with spin a at site i, stands for the nearest 
neighbor pairs of atoms, and t ^ 3 eV [3]. The second 
term accounts for e-e interactions. The expectation val- 
ues of the spin- resolved density = c\^Cicr depend on 
the eigenvectors of %. Thus, a self-consistent field pro- 
cedure is used for solving the problem. The magnitude 
of the on-site Goulomb repulsion parameter U /t ^ 1 was 
estimated on the basis of first-principles calculations and 
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FIG. 1: (color online), (a) Atomic structure of a model of chi- 
ral graphene nanoribbon. The structure shown corresponds 
to 9 — 10.9° chiral GNR characterized by edge repeat vector 
(4,1) and width w — 6. Zigzag and armchair units of the 
edge are indicated, (b) Chirality angles 9 of the considered 
(n + 1,1) and (n -I- l,n) series of chiral GNRs. 
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FIG. 2: (color online). Evolution of the tight-binding band structures (no e-e interactions) of graphene nanoribbons (w = 12) 
upon the change of chirality from zigzag {9 — 0°) to armchair {9 = 30°) via the series of intermediate chiral configurations. The 
scales of the plots account for the varying Brillouin-zone dimensions. The insets show the relation between the one-dimensional 
band structures of zigzag and armchair GNRs and the two-dimensional band structure of ideal graphene. 



experimental data [17]. However, we point out that this 
effective parameter may also incorporate environmental 
factors such as the enhancement of screening due to the 
proximity of dielectric or metallic substrate [13| . For this 
reason we also study the dependence of results on U/t. 

Two parameters determine the structure of smooth 
GNRs: (i) the crystallographic direction of the edge, and 
(ii) the width. In general, the direction of the nanorib- 
bon's edge is defined by some translation vector (n,m) 
of the graphene lattice. For high-symmetry zigzag and 
armchair edges these vectors are (1,0) and (1,1), respec- 
tively. Provided that n > m, the edge of (n,m) nanorib- 
bon is a repeating structure composed of n — m zigzag 
units and m armchair units as illustrated for the par- 
ticular case of a (4, 1) GN R in Fig. Ufa). The length of 
repeat vector a = {n,m) = a^^/v? + nm + m^, where 
ao = 0.246 nm is the lattice constant of graphene. Alter- 
natively, chirality can be described by the chirality angle 

6 = arcsin W | ( , , — j ] . Zigzag and armchair edges 

are characterized by 6* = 0° and 9 = 30°, respectively, 
while for chiral edges 0° < 6 < 30°. We consider chi- 
ral GNRs defined by edge translational vectors {n + 1,1) 
and {n,n + 1) {n> 1). These two series cover the whole 
range of chirality angles 9 and converge io 9 — Q° and 
9 = 30°, respectively, with increasing n [Fig. (rfb)]. The 
configurations of chiral GNRs considered in this study 
cover 4.7° <9 < 25.3°. The width W = y/3waocos9 of 
a GNR is defined by vector {—■w,2'w) pointing along the 
armchair direction as shown in Fig. [Ija). 

We start our discussion by considering the evolution 
of tight-binding band structures of GNRs (neglecting e- 
e interactions, U/t — 0) upon the change of chirality 
from zigzag {9 = 0°) to armchair {9 ~ 30°) via the se- 
ries of intermediate chiral configurations at fixed nanorib- 
bon width (w = 12). High-symmetry zigzag nanoribbon 
[Fig. [21 left panel] exhibits a flat band at the Fermi level 
{E = 0) which spans one-third of the ID Brillouin zone 
(BZ), that is, the corresponding density of edge states 



per edge length per spin p{9 = 0°) ~ l/(3ao). Armchair 
GNRs [Fig. [21 right panel] are either metals or semicon- 
ductors with no electronic states localized at the edges. 
This result can be rationalized within the infinite-width 
picture by considering the 2D band structure of graphene 
projected onto the ID BZ of periodic edge. In the case 
of an armchair GNR both Dirac points {K and K') of 
graphene's band structure are projected onto the F point 
of ID BZ. However, for a zigzag GNR points K and K' 
are projected onto k = 27r/(3a) and k — — 27r/(3a), re- 
spectively, with zero-energy fiat band connecting these 
two points [Fig. [21 insets]. The band structures of chiral 
GNRs in the infinite-width limit can be obtained by con- 
tinuous rotation of the band structure of graphene which 
leads to the known result [l^]: 

pW-^cos(0+^). (2) 

That is, in the limit of large width only armchair GNRs 
show no edge states. The density of edge states is largest 
for zigzag nanoribbons and shows almost linear depen- 
dence on 9. The tight-binding band structures of w = 12 
chiral GNRs [Fig. [21 the x-axis scales correspond to the 
ID BZ dimensions] confirm this picture. Two impor- 
tant comments should be made: (i) in finite-width GNRs 
with edge orientation close to the armchair direction zero- 
energy edge states are partially or even completely sup- 
pressed; (ii) in GNRs with edge orientation close to the 
zigzag direction the flat edge-state band spans whole ID 
BZ and becomes multiple degenerate due to band folding. 

We now discuss the effects of e-e interactions on the 
electronic spectra of chiral GNRs. Figure [SJa) shows the 
band structure and the density of states (DOS) plot for 
a zigzag GNR {w — 12) obtained within the mean-field 
Hubbard model {U/t = 1) compared to the tight-binding 
results {U/t = 0). In the presence of e-e interactions 
the electronic ground state of the zigzag GNR exhibits 
an interesting magnetic ordering: ferromagnetic (FM) 
along the edges and antiferromagnetic (AFM) across the 
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FIG. 3: (color online), (a) Effects of the electron-electron 
interactions on the band structure (left panel) and the density 
of states (right panel) of a zigzag GNR {w = 12). Dashed and 
solid curves correspond to the tight-binding model {U/t = 0) 
and the mean-field Hubbard model {U/t — 1), respectively. 
(b),(c) Respective plots for (2,1) and (6,1) chiral GNRs. 



k = Til a [Fig. late) for d = 7.6° (6,1) GNR]. 

Figured] summarizes the calculated magnetic moments 
and energy splittings for GNRs of various chiralities and 
widths, and different magnitudes oi U/t. Ai U/t = 1 
spin-polarized ground states of few-nm-wide GNRs span 
almost full range of chiralities [Fig. lU^a)]. The mag- 
netic moments per unit length M follow closely the 
dotted curve which corresponds to the complete spin- 
polarization of edge states in the infinite-width limit 
dl]), but appear to be negatively shifted by nearly con- 
stant amounts which are moderately dependent on width. 
Splitting shows a similar dependence on 9 and it is 
independent of w [Fig. HJ^b)]. On the contrary, the mag- 
nitude of A° is largely insensitive to the variations of 
9 in broad ranges of this parameter, but shows a clear 
dependence on w analogous to the case of zigzag GNRs 
[l9|. As the chiral angle 9 approaches 30° magnetic mo- 
ments vanish and the tight-binding band gaps Ag quickly 
rise reaching their maximum values for the correspond- 
ing armchair GNRs [22| . Figure reveal a moderate 
dependence of M{0) on the strength of e-e interactions. 
However, the magnitudes of A'' and A^ are both approx- 
imately proportional to U/t [Fig.HJd)]. 

Finally, we turn our attention to the spatial variation of 
electronic spectra of magnetic graphene edges in relation 



(a) 



nanoribbon [J]. The magnetic moment per edge unit 
length M = 0.313/j,B/ao. Spin-polarization lifts the de- 
generacy of edge states and opens an electronic band gap 
A°. While the tight-binding DOS has only one van Hove 
singularity related to the presence of ID edge states at 
E = the mean-field Hubbard model solution shows two 
pairs of peaks split by A° and A^. We note that Hamil- 
tonian ((TJ respects electron-hole symmetry. Splitting A° 
is related to magnetic correlation between the opposite 
edges while the larger splitting A-^ is due to the FM cor- 
relation along one single edge of most strongly localized 
electronic states at k = ir/a [1, 0, 21 1. Thus, the two 
splittings A° and A^ constitute independent signatures 
of magnetic ordering across the nanoribbon and along its 
edge. The rest of electronic spectrum {\E\ > O.lSi) is 
negligibly affected by the Hubbard term. Van Hove sin- 
gularities at \E\ 0.2t and higher energies correspond to 
the bulk-like states subjected to quantum confinement in 
ID GNRs of the given width. The mean-field Hubbard 
model electronic spectra of chiral GNRs show all the fea- 
tures characteristic of zigzag GNRs provided the ground 
state is spin-polarized. Both splittings, especially A^, are 
reduced in the case of (2,1) GNR {9 = 19.1°) due to the 
smaller magnetic moment M — 0.096^B/ao [Fig. Eljb)]. 
Nanoribbons with chirality close to 6* = 0° show addi- 
tional pairs of van Hove singularities due to splittings of 
the multiple degenerate edge-state bands at A: = and 
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FIG. 4: (color online), (a) Magnetic moment per edge unit 
length AI as a function of chirality angle 9 for three differ- 
ent nanoribbon widths w from the mean-field Hubbard model 
{U/t = 1) calculations. The dotted line shows magnetic mo- 
ments in the limit of infinite width [Eq. ([2])]. (b) Electronic 
band gap A" and maximum energy splitting as a func- 
tion of 9 for different values of w (mean-field Hubbard model, 
U/t = 1). Crosses indicate the tight-binding band gaps Ag. 
The values of (c) M, (d) A° and A^ obtained using the mean- 
field Hubbard model at different values of U/t and u; = 6. 
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FIG. 5: (color online). Variation of the local density of states 
(LDOS) (a) across the (6,1) chiral GNR (x-axis is oriented 
along the edge; y = corresponds to the outermost edge 
atom) and (b) along its edge obtained using the mean-field 
Hubbard model [U/t = 1). (c),(d) Log-linear plots of LDOS 
at the energies indicated in panel (a) across the (6,1) chiral 
GNR and along its edge, respectively. The dashed lines cor- 
respond to the tight-binding LDOS at iJ = 0. The inset in 
panel (d) superimposes the edge structure with the plot. 



to the experimental observations |13|. As a case study 
we investigate the mean-field Hubbard model {U/t — 1) 
local density of states (LDOS) evaluated across the (6,1) 
chiral GNR (w = 12, = 5 nm) [Fig.ISIa)] and along its 
edge [Fig.[5fb)]. We find that both pairs of contributions 
to the total DOS due to the edge states seen in Fig. [3][c) 
are localized at the edge. However, the higher-energy 
LDOS peak aX, E = O.li decays very fast being confined 
within the 1-nm-wide edge region while the lower-energy 
feature dX E — 0.035t penetrates deep into the middle of 
GNR. The series of peaks at > 0.2t repeated in energy 
by « O.lt correspond to the bulk-like graphene states 
subjected to quantum confinement. The total density of 
edge states from the tight-binding calculations shows pro- 
nounced oscillations along the GNR edge [dashed curve 
in Fig. [SJd)]. The oscillation period corresponds to the 
edge periodicity a = 1.6 nm. The E = Q.\t LDOS peak 
from the mean-field Hubbard model follows this trend. 
In contrast, the lower-energy peak {E = 0.035t) shows 
weak variations of LDOS. This behavior was found to be 
generic to all studied chiral and zigzag magnetic GNRs. 

To summarize, our model Hamiltonian study of the 
electronic structure and magnetic properties of chiral 
graphene nanoribbons revealed a number of structure- 
property relations. The described relations can serve 
as unambiguous signatures of edge-state magnetism in 



graphene nanoribbons and provide an avenue towards 
controlling their magnetic and electronic properties. 
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